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With 3 figures in the text 


SUMMARY 


we get, also for the negative moments, the approximate formulae earlier obtained by von 
ardel [2] and Walen [3]. 

' Again starting from Waller’s solution we derive an asymptotic formula, valid for large values 

of a variable x, which is the collision number multiplied by the time elapsed since the entering 
of the neutrons. For small neutron energies we obtain a formula, which can be directly com- 
a ared with those derived by Marshak [4] and von Dardel [2] by means of the method of mo- 
ments. We then see that, although these formulae are rather good approximations in the region, 
where the distribution function is large, they are correct only to the first order if expanded 
in powers of 1/2. = 
os Then in a diagram we give the results of numerical calculations for deuterium and carbon 
_ moderators based on a formula first derived by Svartholm [5]. As Waller [1] has shown, this 
formula follows, in a natural way, from the general solution, and it appears to give, for any 
finite time, the correct asymptotic distribution function for neutrons energies, much smaller 
than the initial energy. In another diagram we compare von Dardel’s [2] distribution function 
and our asymptotic function with Svartholm’s function. Von Dardel’s function is here seen to 
_be a fairly good approximation in the experimentally most interesting region. 
We also calculate the second non-vanishing term of the expansion obtained for the distri- 
bution function by Olsson [6] in the form a power series in time. In order to estimate the 
rate of convergence we then compare this term with the first non-vanishing term, which has 

been calculated by Olsson himself. 

At last we very briefly discuss time distribution functions for different energies. 


Introduction 


At the thime t=0, Q, neutrons with the definite energy HZ, = 410 (vq = velocity, 
-mass=1) enter an infinite homogeneous and_ isotropic moderator consisting 
of free atoms at rest. The neutrons are slowed down by elastic collisions with 
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the atomic nuclei. Capture is neglected. The scattering of the neutrons is as- 
sumed to be isotropic in the center-of-mass system. 
We introduce the following notations: 


M-=mass of moderator atom, measured in units of the neutron mass; 
r=(M—1)/(M+1); 
=(M+1)?/4M; 

du =2 In [((M+1)/(M—1)]= ary; 

v= velocity of a neutron; 

u =In (vp/v") = logarithmic energy variable (lethargy) ; 

l(u) =scattering mean free path of a neutron, assumed to depend on w in 
the followmg way 


L(u)/v=t er" 
where t=1(0)/v, 


and x is a positive constant. (In particular if x= 4, we have J(u) const = 
=1(0)); 
a=vt/l=(t/r)e*"; 
V(u, t)=the average number of neutron collisions per unit w-interval and unit 
time interval; 
N (u, t) = (l/v) ¥ (u, t)=the neutron distribution function. 


Under the assumptions made above we then have Waller’s [1] solution (Wal- 
ler, eq. (24)) 


71 +100 i) +ioco 
Orel PD ee wont 
Yu, )=— | der pq) | dweer-ove “hae fa jcease 
a © c | Fix (1) 
T a ; 2 p(n’) 
T[(l+ +x)/x] 
where = : 
is P((l+y+%— a)/x] Px (7) ea) 
with — 7 (14 oP la duly +1 + 2x) 
P, (n) TL ( tag Phra ; (2b) 


p(n) can also be defined in the following equivalent way (Waller, eq. (12) and (9)): 


= po Cele) - a/n% ’ 
p(n) =e IL @(y + nx) e*'"™, C= Euler’s constant, (3a) 
with G (y) =1—-—*— (1 =e tent 
(helo lena (3) 


G(n) is an analytic function and has one zer 


half plane Re{y} <0. o at »=0 and the others in the 


nO a ad =T(st+1)(re™”)**! if s>0. 
3 ; 


_ 


hus we get by changing the order of the integrations (if s<0, we have to 
a number of integrations by parts in t, before we can change the order of 


- integrations) — 
8 = 3 rai Koes oo 
aac Nf =QrT (¢+ 1) 5-5 | ane pin) | deme 
Os m—ico 0 
eZ apap etn rare p(n) 
Z are: dn pep ee ans | ee Gala (4) 
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If s is a non-negative integer (s=n>0) we now conclude from (4) with the 
help of (8a) that 


ee oo nitioo a 
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If s is a negative integer (s= —m<0), the expression in the right membrum 


of (4) does not exist, but we can calculate its limit value (which exists if 


—u>(m—1)qmu, as we will show later). 
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If u>(m—1)qy, the integration path of the first integral of the aii8 sean 


can be closed to the left, since, according to (3b), we have * 


m-1 j ; 
eT] G(o—kx) =O (eet -Daml) if Re {0} <9. 
K=1 


The integrand has no poles, and consequently the result of the wey i wa is 
zero. Thus from (6) we obtain the result 


ee) Qi tice 


if (—1)""* Cut 1 | D p' (9 —%) x“ 
dt thi" = ; = 
| ¥ (u, ) CG = Thi ae Qn% doe plo TL Ge- kx) 
Q,—100 
(w>(m—I1)qu) (7) 
fatvu,ne 
The time moments GP) (8) 


fat (u, t) 
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can now be calculated from (5) and (7) and we get 
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: vACE n= L Nis ate"), (11) 
f dx (u, x) x” 
Se oe pee (12) 


fdzou, x) 
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_ Using the definition of x and (8) we obtain from these equations 


on ; ‘A ‘ = ; a t?-1 
j 0 ee ae race (13) 
For large u, the equations (10) are applicable, and we put those expressions 
— into (13), getting x 
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appearing in (12) is convergent. Using the logarithmic derivative of (3a) we 
obtain 
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and the right membrum converges under the same conditions as 
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Since, as follows from (3b), 
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In (u’), and consequently also I,,(w’), converge if u’>0, and diverge if u’<0. 
Hence <x ™) exists if u'>0, i.e. if w>mqy but does not emist if u’<0, we. 
if u<mgqy. 
We now see that in order to justify the approximation leading to the sec- 


ond equation (14) we must assume not only large w but also large w’, ice. 
u> Ny. 


q gos ae 


sing the tee sca path to the rahe But the integrand has no poles 
consequently the result is zero. This means that we can use wu as upper — 
t oe deintworal in (15). Thus we get 
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For large values of x the expression 
q ese = (18) 
in (17) can be closely approximated he 
exp {—a[x(u—w) +402 (u—w)? +38 (u—w)]}, (19) 


- 


_ since (19) is a good approximation for (18) if (w—w) is small and for larger 
- (w—w) the expressions (18) and (19) are both so small that we can neglect the 
contributions to the rao from this region. Then we have for large values of x 
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a (1) The method of treating the triple integral (1) which is used here, and the approximation 
_ (19), are due to Professor Waller. 
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(cf. Waller, eq. (53a)). piper Tae 
The integral in the right membrum of (21) can be evaluated according to 
the calculus of residues. If we assume large u, we can neglect all residues except 
the one of the origin, and using (3a) and (3b) we obtain for integral x 
™ 


—xu 2 4 <a , 
pr. Plex) 1st 1 | 
with =" 40) p00) E ea G (kx) (23) 


E=G' (0) 


(cf. Waller, eq. (54)). This is our asymptotic expression for large x and w. 
Now we are going to introduce an approximation for p((z—1)x) into (22) in 
order to get an equation which can be compared with Marshak’s [4] and von 
Dardel’s [2] formulae obtained by the methods of moments. 
From (2b) we conclude that p,((e—1)%)=1 up to any order of 1/2, if x is 
large enough. (!) Inserting this in (2a) we get 


pe—1y9=P (+2) /T (442 (24) 


(*) For light moderator nuclei e °@=[(M-—1) [(M+ ny} is small and then every factor in the 
product of (2b) is very close to unity even for rather small a but for heavy nuclei e “Mal 
and we need very large x-values in order to justify the approximation e @@@**» — 9, 
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and (26) we obi, using, the above formula i re the following 


4 which is correct to the first order in l/s. 
28 , ’ Vu t) r= % get Pa) 27" 
2 sp hace at rah ve 
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with b=4a(e—1) =". (28) 


The right membrum in (27) can now be compared with Marshak’s [4] function 
F )(z) (Marshak, eq. (50), Table I). 

ivi Table 1 we compare our value of 6 with those obtained by Marshak and 
Eon Dardel [2] (von Dardel, Table 1) for M=2, 9 and 15. 


x Table 1 
; M b=(M’-1)'/4M’ ~ Marshak’s b von Dardel’s b 
4 2 0.56 4.86 0.903 
2g ) 19.75 50.15 36.0 
ao 12 35.50 —_ 66.2 (1) 
‘4 15 55.75 128.13 105.1 
be 
A (1) This value has been obtained from von Dardel’s table by means of interpolation. 
3 Marshak’s second approximation F, (x) (Marshak, eq. (54)) cannot be so good, 
since if e*F,(x) is expanded in powers 1/z, which can be made abt successive 
q integrations By. parts, this expansion begins with a*-1 instead of a* 
4 
¢ Numerical results for M=2 and M=12 
ey, 5 -_—- 
a I have made numerical calculations for M=2, «=0(1)14 and M=12, 


= 0(1)14(2)22 using Svartholm’s [5] formula (Svartholm, eq. (6.3)) derived for 


x= and large u, 


~ 22 


15 


Fig. 1. The distribution functions y (x) for M=2 and M= i 


p(x) =e et Stn bE (a), LZ, (x)= Laguerre polynomials 


(cf. Waller, eq. (42)), taking the first 11 terms for M=2 and the first 16 terms 
for M=12. The coefficients were obtained from Svartholm’s eq. (4.9), (4.13) 
with b,=2/&. The L,(x) were taken from tables in an article by Tricomi [7]. 
The result is indicated by the curves of Fig. 1. Our result for M=12 agrees 
with other calculations based on Svartholm’s formula [8]. 


In Fig. 2 we see that for M=2 our asymptotic function obtained from (26) — 
and (11) 


y (x) mee “Too (: ab ib ) 


[shan 162° 


gives a very good approximation to the right of the maximum. But in the z- 
range considered here our asymptotic formula for @ (zx) in the case M =12 shows 
no correlation with the correct formula (cf. note, page 8). 

From Fig. 2 we also see that von Dardel’s function 


7 (x) = 4 gel e 7 blz 
(with A such that [ p(x) dx=1) is a very good approximation for M=2. This 


0 
is the case also for M=12, but I have not drawn the corresponding curve. 
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‘The Seen umes eee in Olsson’ 'S power series 


Olsson, [6] has derived in the case when x=} an pas for (wu, t) in 
t 1e form of a power series in f (Olsson, eq. (32)). 
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4 He has shown Rt <u/qu<m-+1, the coefficients of ¢°, t, t?, ... ’” vanish. 


fe has calculated the first non-vanishing coefficient (Olsson, eq. (35), for which 
he obtained 


= QayrS _ m+Im+6 m+3_, ee 
AE mab pen year i Meer) exp meas ais xy on + (ie + ) 


ages 


ROP 


‘ Aden 


is equal to 


(u'=(m+1)qu—%). (30) 


Further he has discussed the structure of the higher coefficients, which are more 
complicated. According to this discussion the second non-vanishing coefficient 
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_ where M as equal to the sum of residues of the function 
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0+ p(w, x); 1<u/qy<2 and 2<u/qu<3 c ian 
p (x) < is 
After some elementary calculations we get ~ sets 
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The ratio Ky,,2t"*?/(Km4it"*1) is in the case of small energies, i.e. large m, 
nearly equal to 


and we conclude that the second non-vanishing term is small compared with 
the first one, if 
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re : i ate a cae be obtai veal from che) Teac integro- 
equation (Waller, eq. (1) or Olsson, eq. (1); now 5 cae ane 


aa : 
aa ice Be Je | du’ PW, +5(Y 90, (37) 


— af 0 


which has the solution (cf. Olsson, p. 144, first formula) 
E ae y 1 t/t t . 
‘ 3 : ¥(u, )=—e ola ra =| if 0<u/gu< 1. (38) 


From this we get by means of (11) 
x p (u, %) = a(x) 0 (u) +H, (u, x), (39) 


A where f — Yo (a) = = ene (40 a) 


(corresponding to unscattered neutrons), and 
yy (0, 2) =2ae*. (41 b) 


To get thé, form of the curves for small x we use Olsson’s rule above, ac- 
cording to which we have 


p(u,x)~ Km” if m<u/qy<mtl. 


I wish to express my sincere gratitude to Professor Ivar Waller for suggesting the subject 
of this investigation and for zach good advice in the course of the work. 


Tae for Theoretical Phasics, University of Uppsala. 


(?) Numerical calculations are made only for p(x) and p(0, x). The other ~(u, a) (dotted 
in Fig. 3) are only roughly estimated. 
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